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The hyperon-hyperon potentials due to a recent SU(3) Nijmegen soft-core 
potential model are incorporated within a relativistic mean field calculation of 
strange hadronic matter. We find considerably higher binding energy in bulk 
matter compared to several recent calculations which constrain the composi- 
tion of matter. For small strangeness fractions (fs 1), matter is dominated 
by NAE composition and the calculated binding energy closely follows that 
calculated by using the hyperon potentials of our previous calculations. For 
larger strangeness fractions (fs 1), the calculated binding energy increases 
substantially beyond any previous calculation due to a phase transition into 
iVXH dominated matter. We also compare bulk matter calculations with fi- 
nite system calculations, again highlighting the consequences of reducing the 
Coulomb destabilizing effects in finite strange systems. 



I. INTRODUCTION 

Bodmer and Witten independently highlighted the idea that strange quark matter, with 
roughly equal composition of u, d and s quarks leading to a strangeness fraction f s = 
—S/A pa 1 and a charge fraction /q = Z/A ~ 0, might provide the absolutely stable form 
of matter Metastable strange quark matter has been studied by Chin and Kerman 

||. Jaffe and collaborators [f|,0 subsequently charted the various scenarios possible for the 
stability of strange quark matter, from absolute stability down to metastability due to weak 
decays. Finite strange quark systems, so called strangelets, have also been considered [j^3|. 
For a recent review of theoretical studies and experimental searches for strangelets, see Refs. 



Less advertised, perhaps, is the observation made in our previous work |9|,[10[| that 
metastable strange systems with similar properties, i.e. fg ~ 1 and /q ~ 0, might also 
exist in the hadronic basis at moderate values of density, between twice and three times nu- 
clear matter density. These strange systems are made out of nucleons (N), lambda (A) and 
cascade (S) hyperons. The metastability of these strange hadronic systems was established 
by extending relativistic mean field (RMF) calculations from ordinary nuclei (fs = 0) to 
multi-strange nuclei with ^ 0. Although the detailed pattern of metastability, as well as 
the actual values of the binding energy, depend specifically on the partly unknown hyperon 
potentials assumed in dense matter, the predicted phenomenon of metastability turned out 
to be robust in these calculations [ TO), 1 1 ] . 
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Quite recently, Stoks and Lee |12| have challenged the generality of the above results for 
strange hadronic systems. These authors constructed G matrices for coupled baryon-baryon 
channels, using an SU(3) extension |Tj| of the Nijmegen soft-core NSC97 potentials [H from 
the S = 0, — 1 sector (to which data these potentials have been fitted) into the unexplored 
S = —2, —3, —4 sector. These G matrices were then employed within a Brueckner-Hartree- 
Fock (BHF) calculation of strange hadronic matter (SHM) in bulk. The results showed that 
NAE systems are only loosely bound, and that charge-neutral strangeness-rich hadronic 
systems are unlikely to exist in nature in metastable form, in stark contrast to our earlier 
findings P,|T0|. 

This vast difference in the predictions for the metastability and binding of SHM between 
following a BHF methodology, which uses an SU(3) extrapolated form of the NSC97 baryon- 
baryon potentials, and following a RMF methodology, which is based on mean fields designed 
to mimic the consequences of the Nijmegen hard-core potential model D [15], has prompted 
us to investigate possible origins of it. In this work we present calculational evidence for the 
incompleteness of the procedure applied by Stoks and Lee 
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We do so by reproducing 
qualitatively their results for the instability and weak binding of iVAH matter in bulk, 
within a constrained RMF calculation in which the mean fields are now designed to mimic 
the consequences of the NSC97 model used by Stoks and Lee. The constraints imposed by 
us, as a check, are identical with those imposed by these authors for the composition of SHM 
(see Fig. 4 of Ref. |12|)- We argue that this is not the right way to identify minimum-energy 



equilibrium configurations for SHM. Indeed, doing the unconstrained RMF calculation with 
the same NSC97-inspired mean fields, we find qualitatively good agreement, for f$ < 1 
where the bulk matter is iVAH dominated, between these new results and our old results 
in model 2 [10|. For fg ^ 1, the new unconstrained calculation results in considerably 
higher binding energies than ever calculated for SHM, due to a phase transition into A^EH 
dominated matter. 

The paper is organized as follows. In section II we describe the methodology of finding 
equilibrium configurations within the RMF formalism, and the input mean fields entering 
the new RMF calculations. Section III includes the results of these new calculations for 
bulk SHM, as well as for finite multi-strange systems for which BHF calculations have not 
been done to date. The role of the Coulomb interaction in stabilizing charge-neutral strange 
systems is highlighted. Our results are summarized and discussed in section IV, where we 
also comment on the applicability of the SU(3)-extended NSC97 potential. 



II. METHODOLOGY AND INPUT 

We adopt the Relativistic Mean Field (RMF) Model to describe strange hadronic matter 
(SHM) in bulk and for finite systems of nucleons and hyperons. The model is an effective 
model where the parameters are adjusted to the known properties of nuclei and hypernuclei. 
We include in our extended RMF model all the l/2 + baryons of the lowest SU(3) flavor octet, 
as well as hidden-strangeness meson exchange to allow for possibly strong hyperon-hyperon 
(YY) interactions. Here we use model 1 and model 2 of Ref. [|§. The basic ingredients of 
these models are the octet baryons matrix B, the matrices V^ 8 ** and Vjp of the vector meson 
octet and singlet, respectively, and the two scalar mesons a and cr*. In addition, a Coulomb 
term is included in finite system calculations. The Lagrangian is given as 
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C = TrB(i^dfj, - g aB o - g a *BO* - m B )B 



^3 C 4 L { a *au * 2 * 2 

V .t- , 3 4^ ~ 2 i ^ '"'.^ 

(aTr^ [l/ (8) , 5] + (1 - a) TrBY {v^ (8) , B}) - g^TrB^B ■ TrV^ 

4* 



--W;^ + -Trmy;^ + 7 dK^) 2 . (1) 



Here, both scalar fields are treated as singlets. The octet vector fields can be coupled in 
two ways, either antisymmetric (F-type with a = 1) or symmetric (D-type with a = 0). 
In the mean-field approximation, only the u, p, and <fi vector mesons remain operative. 
Their coupling constants to the baryon fields can be related by SU(3) symmetry [TjJ. By 



assuming ideal mixing of the vector mesons (i.e. the is a purely ss state), pure F-type 
coupling (a = 1), and that the nucleon does not couple to the 4>, one recovers the SU(6) 
relations of the simple quark model 

1 1 1 

-QujN — ~9ujK — -QujY, — 3wS 

9pN = 2#p£ = 9pE , 9 P A = 

2<?0A = 2#0s = g<f,3 = — —g^N , g<j,N = . (2) 

Here, the constraint 3g p ^ = g^N is relaxed to allow the isovector coupling constant to be 
fixed by the isospin dependence of nuclear binding energies. The quark model is not used 
for the scalar coupling constants, rather they are determined by adjusting to nuclear and 
hypernuclear properties. Self-interaction terms for the scalar field a and the vector field u 



are also included in the model. In the following, we use the parameter set TM1 of Ref. fl7 
where the parameters were taken from a fit to properties of spherical nuclei. The remaining 
scalar a coupling constants for the hyperons are chosen to give reasonable hyperon potentials 
in saturated nuclear matter: 

U^\ Po ) = -30 MeV , U^\p ) = +30 MeV , and ujpfa) = -18 MeV . (3) 

Note, that these relativistic potentials are (10 — 20)% stronger than the corresponding non- 
relativistic values. The values for E and H hyperons differ from the previous choice of Refs. 
PJT0[|, reflecting recent developments in hypernuclear physics which are briefly recorded 



below. 

For the E nuclear interaction, the most updated analysis of E~ atomic data indicates 
a repulsive isoscalar potential in the interior of nuclei [pi which is compatible with the 
absence of bound-state or continuum peaks in a recent search for E hypernuclei |TJ . In 
fact, the only E hypernuclear bound state found so far is |He pO| , ^T| , where the binding 



results from the strong isovector component of the E nuclear interaction. These statements 
are supported by several recent calculations p2] , |23[| . The precise magnitude of the depth of 



the E nuclear potential is of little importance to our investigations. It turns out that the 
E hyperon will not appear anyway in the bulk matter calculation, or in the finite system 
calculations within models 1 and 2, unless its hadronic interactions are exceptionally strong, 
so as to block the release of about 75 MeV in the free-space EB — ► AB strong-interaction 
conversion. 
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For the H nuclear interaction, measurements of the final-state interaction of S hyperons 



produced in the (K ,K + ) reaction on 12 C in experiments E224 at KEK [24] and E885 at 



the AGS p5[ indicate a nonrelativistic potential Ug^ of about —16 and —14 MeV or less, 

respectively. Below we will actually vary the value for to check its effect on the binding 
energy of SHM. 

The hyperon (Y) potentials U Y in hyperon (Y r ) matter, in the absence of direct ex- 
perimental data, depend to a large extent on the assumptions made on the underlying YY 
interactions. In model 1, which does not use a* and exchanges, the potentials Uy are 
rather weak, less than 10 MeV deep. The exchange of these hidden- strangeness mesons is 
included in model 2, where the a* coupling to hyperons is adjusted so that the potential of 
a single hyperon, embedded in a bath of H matter at nuclear saturation density po, becomes 

t4 S) (p )=t/P(po) = -40MeV , (4) 

in accordance with the attractive YY interactions of the Nijmegen potential model D ||10|| . 
The resulting U^ is about —20 MeV, considerably more attractive than in model 1. Indeed 
the few double A hypernuclear events observed so far in emulsion require a relatively strong 
AA attractive interaction p6| , which lends support to model 2 over model 1, but the actual 
situation for the other, unknown, YY' channels could prove more complex than allowed 
for by either model. All that may be said at present is that, as far as the AA interaction 
strength is concerned, model 2 is a more realistic one than model 1. 

Since there appears some confusion in the recent literature [0,^] regarding how to 
calculate self consistently the properties of SHM in bulk, we will ponder on the thermody- 
namically consistent methodology in more detail. As we will demonstrate in the following, 
a major property of SHM within the SU(3)-extended NSC97 model might have been over- 
looked in these works. Here we focus on the thermodynamically correct treatment in the 
RMF approximation. The extension to BHF calculations is then straightforward. Very 
recently, thermodynamically consistent BHF calculations of /3-stable strange matter in neu- 
tron stars have been performed by Baldo et al. ||28| , using the NSC89 model [2S] for the YN 



interactions, and by Vidana et al. f30|, using the SU(3)-extended NSC97 model |l3j for the 
YN and YY interactions. 

In general, we can describe the system by the grand-canonical thermodynamic potential 
Q, which depends on the temperature T, the volume V, and the independently conserved 
chemical potentials fj, a . At T = 0, the pressure is given by: 

P(ji a )=-n(ji, a ,T = 0)/V . (5) 

For SHM in bulk, since the isospin dependence is usually suppressed, there are just two 
conserved charges in bulk which are the baryon number B and the strangeness number S. 
The chemical potentials of the individual baryons can be related to the corresponding baryon 
chemical potential hb and strangeness chemical potential \is by 

Hi = Bi ■ hb + Si ■ [i s • (6) 

This ensures that the system is in chemical equilibrium or, in other words, that the 
strangeness and baryon numbers are conserved in all possible strong-interaction reactions in 
the medium, such as 
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E + N ^ A + N A + A^S + iV A + S^S + S ... (7) 

The Hugenholtz - van- Hove theorem relates the Fermi energy of each baryon to its chemical 
potential in equilibrated matter 



Hi = Epj = yj k 2 Fi + m* + g^uo + g pi p + g&fa . (8) 

Here we used the energy-momentum relation of the mean-field approximation with the ef- 
fective mass m* = mi + g ai cr + g a *i(J* for the baryon species i. Note, that these potentials 
depend on the overall composition of the matter, requiring thus a self-consistent calculation. 
Equation (|8|) can be easily solved to calculate the Fermi momenta kp^ and hence the number 
density of baryon i for given chemical potentials, 

k 3 ■ 

where ji is the spin-isospin degeneracy factor. If the solution results in an imaginary Fermi 
momentum, the particle is not present in the system and the corresponding density is set to 
zero. In Brueckner theory, one has to solve for an equation of the form 

k 2 

^ = Ei{k F>i ) = m i + ^ + toUi{k F ,i) (10) 

since the potential is now momentum dependent. It is apparent from the above procedure 
that no baryon species may be ignored a priori, but one has to check for their appearance 
by calculating the corresponding Fermi momentum. Therefore, any calculation of baryonic 
matter with nucleons and A's alone is bound to violate the condition of chemical equilibrium 
Eq. (0), since E and particularly 5 hyperons are likely to appear in strange hadronic systems 



||1 1|| . In that sense, multi-A matter calculations as performed in |31|-p^1 are incomplete. 
Furthermore, the Fermi momenta of different baryon species cannot be set equal to each 
other, since this again violates the condition of chemical equilibrium. If one arbitrarily sets 
certain baryon fractions to be equal to each other, as e.g. done in ||12|| , the resulting system 



is not in its energetically favored global minimum and the computed binding energies will 
be underestimated. In addition, the overall pressure of the system will be too low, resulting 
in a too soft equation of state; maximum masses of neutron stars computed in this way will 
be underestimated. 

The pressure and the energy density in the RMF model in bulk are given by 

P = ~\mla 2 - -a 3 - - A a A - -m^a* 2 + im> 2 + -d^ + \m\<& 
r k F,i „ k 2 



i=B ' jZ (2tt) 3 7o ' ^k 2 + mf 

loo b o C a Xo 2 1 22 ^ 7 A 2/2 

e = -ran H — o -\ — a H — m„,<7 H — m, ,u n -\ — au n -\ — m.cbn. 
2 3 4 2 ° 2 w0 4°2^ 

respectively. The binding energy per baryon is then obtained by subtracting the properly 
weighted combination of the rest masses from the energy density of the system 
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E/A = — (e - p N ■ m N - p A • m A - p E • m s - p H • wis) 
Pb 



(12) 



The meson fields are determined by their equations of motion (see e.g. [38]). The particle 
densities are calculated using the thermodynamically consistent formalism as outlined above. 

For most purposes, one wishes to translate the dependence on the two chemical potentials 
into their corresponding baryon and strangeness number density. This is done using the 
expressions for overall baryon and strangeness number conservation 



Pb = J2B { - p V7i = p N + p A + p E + p s 

i 

PS = Y;(~ S i) ' Pv,% = PA + Ps + 2p s 



(13a) 
(13b) 



where pv,i is the vector density of baryon species i. 
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FIG. 1. Binding energy per baryon of SHM in models 1 and 2, with different choices of the H 
potential in nuclear matter. 



A useful measure of the strangeness contents of the system is given by the strangeness 
fraction: 



f Ps 
Is — — 



Pa + Ps + 2p s 



Pb Pn + Pa + Ps + Ph 



(14) 



III. RESULTS 



Plots of calculated binding energy of SHM per baryon as function of the strangeness 
fraction fg are shown in Fig. [I], for different choices of the H nuclear potential denoted by 
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Us, in models 1 and 2 |||TIJ. For fs = 0, there are only nucleons in the system and one 
gets the standard equation of state of nuclear matter as function of baryon density. The 
equilibrium density of nuclear matter is determined by minimizing the binding energy with 
respect to the baryon density. The resulting minimum value of binding energy per nucleon 
is shown then at fs = in the plots of Fig. Next, we increase the strangeness fraction 
from zero on, and the system of equations adjusts itself at each fixed value of fs to find the 
corresponding baryon densities ensuring chemical equilibrium (Eq. (|5|)). The minimum value 
of the binding energy per baryon as function of baryon density at each fixed strangeness 
fraction is then plotted in Fig. [TJ for the corresponding value of fs- In this way, one gets the 
binding energy of SHM as function of the strangeness fraction. It turns out that E hyperons 
do not appear at any value of fs in both models 1 and 2. To display the dependence on the 
S nuclear potential we chose three different values, Us = —10, —18, —28 MeV. The variation 
in the plots of model 1 is quite pronounced. For Us = —28 MeV, the minimum is at a 
finite value fs = 0.6, with a binding energy per baryon of —17 A MeV. For shallower H 
potentials, this minimum disappears and slightly strange matter with fs ~ 0.1 is the most 
strongly bound configuration. On the other hand, in model 2, varying Us does not lead to 
drastic changes. The minimum in the binding energy per baryon for Us = —28 MeV, at 
f s = 1.3 with E/A = -24.6 MeV, is shifted to E/A = -21.5 MeV for U B = -18 MeV and 
to E/A = —19.6 MeV at a slightly higher value fs = 1.4 for Us = —10 MeV. The reason is 
that in model 2 the minimum is generated by the YY interactions which have been adjusted 
according to Eq. (|4]), so that the binding energy curves in model 2 are not as much affected 
by changing Us as compared to the effect of this change in model 1. Note that the constraint 
(|j) ensures that pure S matter (fs = 2) has the same binding energy, E/A = —8.9 MeV, in 
all three cases. Pure H matter is always unbound in model 1 due to the missing attraction 
in the YY channels. 

Substantial departures from the universality (Eq. @) assumed in Refs. [§,|H| for the YY 
interactions occur in the most recent SU(3)-extension of the Nijmegen soft-core potential 



model NSC97 |13| . In particular, the EE and SS interactions are predicted to be highly 
attractive in some channels, leading to bound states. We wish to examine the consequences of 
this model in our RMF calculation of SHM. The YY interactions of Ref. |13j are implemented 
in our calculation by adjusting the coupling constants of the a* meson field to reproduce 
qualitatively the hyperon binding energy curves shown in Fig. 2 of Ref. for set NSC97f. 
All the other coupling constants are held fixed, so that we still get the hyperon potentials 
of Eq. (^) in nuclear matter. The resulting binding energy curves, of each baryon species 
j in its own matter are depicted in Fig. |2| as function of density. For nucleons, we 
again use the parameterization TM1 so as to get the correct binding energy at the correct 
saturation density p . Note that NSC97f does not reproduce the correct nuclear matter 
saturation point, but gives a too shallow minimum at a too high density (see Fig. 2 of Ref. 
||12||). No binding occurs for A hyperons, and B^ reaches +20 MeV already at rather low 
density, p = 0.1 fm -3 . This strong repulsive "potential" is due to the very weak underlying 
AA interaction in the extended NSC97f model which is incompatible with the fairly strong 
A A attraction necessary to explain the observed double A hypernuclear events (see pB,3D 
and references therein). On the other hand, E matter is deeply bound, by —33 MeV per 
baryon at p = 0.58 fm~ 3 which is twice as deep as ordinary nuclear matter, and S matter 
has a binding energy of —23 MeV per baryon at p = 0.39 fm -3 . It is clear from Fig. that a 



7 




Density p (fm 3 ) 



FIG. 2. Binding energy per nucleon (N) in nucleon matter, compared to the binding energy per 
hyperon (A, E, 3) in its own hyperonic matter. The hyperonic parameters were chosen to reproduce 



the binding energy minima of Fig. 2 in Ref. [12]. 



mixture of E and H matter must be very deeply bound too, unless there is an overwhelmingly 
repulsive interaction between E and 5 hyperons. Actually, the interaction between E and 
H hyperons is the most attractive one in the extended NSC97f model, giving rise to the 
deepest bound ES dibaryon state f|X3f1 . Independently, increasing the number of degrees of 



freedom will also result in a more deeply bound state. This is the case, for example, when 
going from unbound neutron matter (7 = 2) to bound nuclear matter (7 = 4). Therefore, 
one expects that ES matter is in fact more deeply bound than E or S matter alone. In the 
following, we will denote the parameterization responsible for the curves of Fig. |2| as model 
N. 

Fig. ^ shows the binding energy of SHM per baryon in model N as function of strangeness 
fraction. For comparison, the curves for model 1 and 2 from Fig. [T] are also plotted. We 
performed two different calculations for model N: one where the hyperon fractions \% = Pi I ' Pb 
are held equal by hand (xa = Xs = Xs) as done in Ref. fll2| , and the self-consistent one 



where the hyperon fractions are determined so as to ensure chemical equilibrium (denoted 
as "equil." in the figure). It is evident from Fig. || that the self-consistent treatment gives 
a substantially lower energy, since it is the unconstrained minimum-energy solution. The 
disagreement between the results of the two calculations increases with fg. At fg = 1 the 
difference between the two curves amounts to nearly 10 MeV. The curve for model N follows 
closely the one for model 2 up to a strangeness fraction of fg = 0.95. At larger values of 
fg a deep minimum develops due to the highly attractive interaction between the E and S 
hyperons. Note that an equal mixture of E's and S's gives fg = 1.5 and the minimum of 
the curve is close to that point, i.e. E/A = —79 MeV at fg = 1.45. For larger values of fg, 
the curve denoted by N rises again, ending up at fg — 2 with the same binding energy per 
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FIG. 3. Comparison of the binding energy of SHM per baryon in models 1 (dash), 2 (dash-dot) 
and N (solid). The upper solid line shows the result for the constrained case of equal hyperon 
fractions (xa = Xs = Xh)> the lower one shows the curve for the correct, unconstrained equilibrium 
calculation. 



hyperon as that shown in Fig. fj for pure S matter. The deep minimum around fs = 1.5 
results from the deep binding of £ plus H matter which is stronger than seen in Fig. ||] for 
matter composed of either one of these species separately. This deep minimum structure can 
only be reached when the hyperon fractions are allowed to adjust self consistently. Fixing 
the hyperon fractions, as done in Ref. W2\, will always give a curve which is higher in 



energy, risking the loss of some important features of the model. In fact, the curve for the 
constrained calculation JTJ! ends up at fs — 4/3 due to the particular constraint applied. 

The deep minimum seen in Fig. [3] emerges due to a second minimum in the corresponding 
equation of state at high strangeness fraction, connected with a first order phase transition 
from matter consisting of iVAS baryons to iVEH baryonic matter. This transition is visu- 
alized in Fig. |] where the binding energy is drawn versus the baryon density for several 
representative fixed values of fs- For fs = 0.8, there is a global minimum at a baryon 
density of ps = 0.27 fm~ 3 . A shallow local minimum is seen at larger baryon density at 
Pb = 0.72 fm~ 3 . Increasing the strangeness fraction to fs = 0.9 lowers substantially the 
local minimum by about 20 MeV, whereas the global minimum barely changes. At fs = 1.0 
this trend is amplified and the relationship between the two minima is reversed, as the min- 
imum at higher baryon density becomes energetically lower than the one at lower baryon 
density. The system will then undergo a transition from the low density state to the high 
density state. Due to the barrier between the two minima, it is a first-order phase transition 
from one minimum to the other. 

Fig. [5] demonstrates explicitly that the phase transition involves transformation from 
NAS dominated matter to A^ES dominated matter, by showing the calculated composition 
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FIG. 4. Transition from NAE matter to iVXH matter in model N. A second minimum appears at 
higher density for a strangeness fraction of fs = 0.8, becoming more stable for higher strangeness 
fraction (fs = 1). 




FIG. 5. Composition of SHM in model N versus the strangeness fraction. S hyperons appear 
around fs = l. 
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of SHM for model N as function of the strangeness fraction fg. The particle fractions Xi 
for each baryon species change as function of fs- At fs = 0, one has pure nuclear matter, 
whereas at fs = 2 one has pure H matter. In between, matter is composed of baryons 
as dictated by chemical equilibrium. A change in the particle fraction may occur quite 
drastically when new particles appear, or existing ones disappear in the medium. A sudden 
change in the composition is seen in Fig. for fs = 0.2 when E's emerge in the medium, 
or at fs = 1.45 when nucleons disappear. The situation at fs = 0.95 is a special one, as 
S's appear in the medium, marking the first-order phase transition observed in the previous 
figure. The baryon composition alters completely at that point, from NE baryons plus a 
rapidly vanishing fraction of A's into SH hyperons plus a decreasing fraction of nucleons. 
At the minimum of the binding energy curve in Fig. |||, the matter is composed mainly of 
E's and S's with a very small admixture of nucleons. 




FIG. 6. Binding energy of multi-strange finite systems built on a Sn nuclear core in models 1 
(circles) and 2 (squares, stars), with and without Coulomb effects. The curves for bulk SHM (solid 
lines) are also shown. 



Last but not least, we discuss finite systems of SHM for which the Coulomb interaction 
plays a significant role. We also performed calculations switching off the Coulomb interaction 
in order to study separately the effects due to the possibly strong YY interactions. Yet, 
we will not pursue for finite systems the implications of the deep minimum found in model 
N, but rather stick in the following to the more conservative models 1 and 2. Our choice 



for the H potential is £/= = —18 MeV, in accordance with the recent observations p4 , 25 



Following the procedure outlined in P,|T0"|| we start from a "normal" nucleus (here 100 Sn, 
for example) and add A hyperons to the system. As soon as the strong decay reactions 
— > AA or riEP — > AA are Pauli-blocked, we start adding 5 hyperons. For each given 
system of nucleons and hyperons, we check that the two reactions are blocked also in reverse, 
so that the whole multi-strange nucleus is metastable, decaying only via weak interactions. 
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These calculations also include the effects of the p meson field, in order to properly account 
for symmetry-energy contributions to the binding energy. 

Our results are summarized in Fig. |6|. The solid lines are the binding energy curves for 
SHM in bulk for model 1 (upper curve) and model 2 (lower curve), taken from Fig. |l| for the 
value Us = —18 MeV. The filled symbols denote systems where the Coulomb interaction 
is included, the open symbols stand for the case where the Coulomb interaction has been 
switched off. For model 1, shown in circles, SHM is slightly more bound than the core nucleus 
100 Sn. For the highest strangeness fraction we found, fg = 0.375, the system is bound by 
E/A — —9.5 MeV, compared to E/A — —8.3 MeV for 100 Sn. When the Coulomb interaction 
is turned off, the curve shifts down by several MeV per baryon. The chargeless core nucleus 
100 Sn is then bound by E/A = —12.2 MeV, whereas the strangeness richest object has now 
a slightly lower binding energy per baryon of E/A = —11.9 MeV. Hence, the main effect 
for the increased binding energy of SHM in model 1 actually comes from Coulomb effects. 
The negatively charged H~ hyperons neutralize the positively charged protons, making SHM 
more bound than ordinary nuclei. Still, the binding energies found are well above the curve 
for bulk matter due to finite size effects, such as surface tension. 

The situation is different in model 2, for which results are shown by squares in Fig. 
|5|. Due to the attractive YY interactions in this model, the binding energy per baryon 
increases substantially to a value of E/A = —19 MeV at fg = 1.3, which is even deeper 
than the binding energy of nuclear matter in bulk. This high value of binding is obtained 
irrespective of whether or not the Coulomb interaction is included, since for such high values 
of strangeness fraction the total charge fraction of the system is close to zero. Obviously, 
in model 2 the tremendously increased binding energy of SHM originates mostly from the 
attractive YY interactions, and only to a minor extent from the reduced Coulomb repulsion. 
Note that the binding energy for the deepest lying systems in model 2 is close to the values 
for SHM in bulk matter which are shown by the lower solid line. These systems are quite 
heavy, with mass numbers of about A = 400 and higher, so that finite size effects become 
quite small. 

In addition, we also plotted the binding energies of purely hyperonic systems which 
consist of AH°H~ hyperons solely and, therefore, do not need to be Pauli blocked in order 
to keep them metastable P,|l0"||. Since these systems can decay only via weak interaction, 
arbitrary numbers for the three different hyperon species are allowed. We find that purely 
hyperonic objects are bound up to E/A = —12 MeV per hyperon. The binding energies, 
denoted by stars in the figure, follow the trend of the bulk calculation curve (solid line). 



IV. SUMMARY AND CONCLUSIONS 

In the present work we have calculated the minimum-energy equilibrium composition 
of bulk SHM made out of the SU(3) octet baryons N, A, £ and S, over the entire range 
of strangeness fraction < fg < 2, for meson fields which generate, within the RMF 
model, qualitatively similar baryon potentials to those generated from the SU(3) extension 



of the NSC97 potential model using the BHF approximation |ji2| . Our main results 
are displayed in Fig. [3] which shows that SHM is comfortably metastable in this model N 
for any allowed value of fg > 0. The NAE composition and the binding energy calculated 
for equilibrium configurations with fg < 1 resemble those of model 2 in our earlier work 
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PJlO| . The use of models 2 and N [f3~Q| ? p~3[] implies sizable Y Y attractive interactions which 
differ, however, in detail between the two models. The model of Ref. yields particularly 
attractive EE, EE and EH interactions, but vanishingly weak AA and A/H interactions. 
We remark that this extent of weakness is in fact ruled out by the little information one 
has from AA hypernuclei [f26] and from H-nucleus interactions [p4| , p5| . On the other hand, 
model 2 of Ref. ]l(J accounts more realistically for the attractive AA and A/E interactions, 
but ignores altogether E hyperons which require exceptionally strong binding in order to 
overcome the strong-interaction Ei? — > AB conversion which in free space releases about 
75 MeV. Yet, all these differences between the two models regarding the relative size of 
interactions within A^AH dominated matter hardly matter when it comes to establishing 
the stability and binding pattern of this multi-strange matter. In this sense, SHM is a 
robust phenomenon. The metastability of SHM has also been recently confirmed within the 
modified Quark-Meson Coupling Model J|0| . 

The difference between models 2 and N clearly shows up for fg > 1, where E's replace 
A's in model N due to their exceptionally strong attraction to E and H hyperons. Figs. [3|, f|, 
|5] of the present work give evidence for a phase transition, from A^AE dominated matter for 
fs 1 to A^ES dominated matter for fs <^ 1, with binding energies per baryon reaching as 
much as 80 MeV. This effect has gone unnoticed in previous works which by constraining 
the composition of matter in bulk did not allow for the most general minimum-energy 
equilibrium configurations. In contrast, our model 2 produces a much smoother pattern of 
binding over the entire range of fs, with a gain of only approximately 5 MeV per baryon 
(at fs ~ 1.3) for the bulk matter calculation. However, for finite multi-strange systems the 
gain can be considerably bigger, due to getting rid of most of the Coulomb repulsion for 
such approximately charge-neutral systems, amounting to almost 11 MeV per baryon for 
the examples of Fig. || 

We checked also for the critical strength of the H-nuclear potential below which finite 
systems of SHM would consist only of nucleons and A hyperons in model 1. Of course, the 
critical value for Us depends on the size of the system. For a nuclear core of 16 with 8 A's 
filling up the s- and p-shells, H's cannot be added to the system for a potential shallower 
than m = -13 MeV. In the case of 56 Ni, this critical value shifts to £7| = -10 MeV. For 
the 100 Sn nuclear core used to demonstrate finite systems of SHM in our present calculation 
(see Fig. |]), we find a critical strength JJ% = —7 MeV for which the SA" — > AA strong- 
interaction conversion is barely Pauli-blocked. However, Us needs to become more attractive 
than the above critical values demonstrate, in order that the corresponding multi-strange 
finite systems also remain particle stable. We remind the reader that the value Us = —18 
MeV used in the figure was designed to agree with the present phenomenological estimates 



2J,|25[. In bulk matter, H-nucleus potentials as repulsive as Us = +40 MeV still admit 
bound H's just before SHM gets unbound at fs = 0.7. The reason for this behavior is that 
the constraint fs = 0.7 introduces H's in matter even though their interaction is repulsive. 
The Fermi momentum of the A's become sufficiently high so that it pays to create some 
seemingly unfavorable H's in order to lower the A Fermi momentum. Note that SHM in 
model 2 always contains H's, irrespective of the nature of Us, by virtue of the attractive 
underlying YY interactions. 

While it is true that the RMF model 2 is a schematic model and is linked only indirectly to 
the underlying baryon-baryon interactions, it is nevertheless constrained by A and H nuclear 
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phenomenology, and by the few AA hypernuclear species reported to date. The extrapolation 
to YY channels which underlie the hyperon potentials in hyperon matter is more conservative 
in this model than in model N inspired by the SU(3) extension of the NSC97 potential model 

13| . We emphasize that although the NSC97 model fl4j has been tuned up to reproduce 
certain characteristics of A hypernuclei, particularly its version NSC97f, the predictions 
elsewhere of these models appear invariably ruled out by whatever experimental hints one 
has to date. In addition to the exceedingly weak AA and iVS interactions already mentioned 
above, the NSC97 model overbinds A hyperons in nuclear matter (U^ ~ —38 MeV) and 
gives rise to quite attractive £ nuclear potential ~ —20 MeV) in BHF calculations 

27J, whereas the phenomenology of S~ atoms [|R| and 'hypernuclei' |22|] indicates a much 



weaker, if not a repulsive, £ nuclear potential. Furthermore, the NSC97 model gives rise to a 
sizable £ nuclear symmetry energy which is opposite in sign [f27| to that of the earlier NSC89 
model 



29 



and, more importantly, to that established phenomenologically |18| , p2| , |23l . We 
therefore suggest that the consequences of the NSC97 model in dense SHM, as exemplified 
here, should be taken with a grain of salt. More dedicated work is required to amend the 
pitfalls of this model (for a recent discussion in this direction see Ref. |41[|). 
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